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Summary. I briefly review the advancements in the construction of the graviton propagator 
in the context of LQG and Spinfoam Models. In particular the problems of the Barrett-Crane 
vertex in giving the correct long-distance limit and the introduction of the new corrected 
models. This kind of calculation applied to an alternative vertex with given asymptotic can 
give the correct propagator and is then able to help selecting spinfoam models. In particular 
the study of the asymptotic properties of the new models shows the predicted behavior able 
to overcome the BC difficulties and to give the correct propagator 

Loop Quantum Gravity (LQG) [I] is a promising candidate for a theory of Quantum 
Gravity based on a canonical quantization of General Relativity (GR). Its dynamics is defined 
through the action of the Hamiltonian constraint [2] but until now, the implementation of 
this constraint has encountered many difficulties that have led to path integral formulations 
of the same theory called Spinfoam Models (SM) pQ. The most promising and extensive SM 
studied in the literature was the Barrett-Crane (BC) model |3] . The main open problem is the 
contact with the low energy world namely on a way to recover from LQG and SM low energy 
predictions described by GR and its linearized quantization. Remarkably in the last few year 
a technique for computing n-point functions in this background-independent context has 
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been introduced in [3] \E\ and developed [5] (see also [7] and references therein). One of the 
main achievement of this technique has been the fact that this kind of calculation has ruled 
out the BC model [8] and has lead to the introduction of a new class of spinfoam models. 
In fact the difficulties found in the construction of the propagator with the BC dynamics 
have focused the attention on the intertwiner dependence of these models (associated to the 
angles between simplices) revealing a strong imposition of the simplicity constraints in the 
quantization procedure that has the effect of freezing some degrees of freedom. This problem 
has been corrected in a new model now called EPR [9] that has the remarkable feature of 
having boundary states compatible with the boundary states of LQG. This vertex can be 
defined for general values of the Immirzi parameter 7 [TO] and extended to the Lorentzian 
case [11] . A similar vertex now called FK |12] has also been derived using the coherent 
states techniques introduced by Livine and Speziale [13]. For 7 < 1 the two techniques yield 
exactly the same theory. These models can now be tested with the propagator technique and 
another complementary based on the wave packet propagation |14l 1151 116] to understand 
their semiclassical limit 17, 18, 19, 20 and to construct Feynman rules for quantum gravity. 

1 LQG propagator: BC model and the New models 

Given a regular 4-simplex with two boundary tetrahedra n and m centered at the points 
x and y we can define the projection of the euclidean graviton propagator G ,lvf "* (x,y) as 
Gtt;%(L) = G M " p<T (a;,y)(nS I i) ) M (n^ ) )„(n^ ) )p(n^ ) ) (T , where the latin indexes label the five 
tetrahedra bounding the 4-simplex and n^' is the normal one-form to the triangle bounding 
the tetrahedra m and k, in the hyperplane defined by m, and L is the euclidean distance 
between x and y. G^;„(L) can be computed [4] in a background independent context as 

G<ft' m = (W\ (£« ■ Eg> ~ ■ n«»)(^ ■ £<» - ■ r&>) |# q ). (1) 

for an appropriate q. We refer to UJ as the LQG graviton propagator. The computation of 
this object is based on three ingredients: the boundary functional {W\, the triad operator E„ 
at the point n, contracted with n„ and the boundary state \&q), picked on a given classical 
boundary (intrinsic and extrinsic) geometry q. The diagonal components Gq 1 ^^ have been 
computed 5 using the BC model to codify the dynamics, the double grasping operators of 
LQG acting in the same directions (the LQG area operators) and a vacuum state given by a 
gaussian superposition of spinnetwork states. In [5] was shown that at large distance, these 
components of UJ agree with the ones of the conventional graviton propagator. 

The construction of the non diagonal terms [8] [21] with the same tecnique, required a 
change in each of the three ingredients, because the graviton operators ■E^ > (associated 
with the 3d dihedral angles between the triangles in the tetrahedra n) called into play the 
dependence of the spinnetworks from the intertwiners and in turns, the dependence of the 
boundary state and the vertex from these variables. In particular the BC dynamics used to 
compute the diagonal terms has a trivial intertwiner dependence that appeared insufficient 
to deal with the non diagonal terms. To see how this works it is enough to calculate (jTJ) to 
first order in the GFT [JJ [55] expansion, and in the limit in which the boundary surface is 
large. Eq. UJ receive the leading contribution for W with support only on spin networks with 
graph dual to a 4-simplex. If j and i are, respectively, the ten spins and the five intertwiners 
that color this graph, then in this approximation UJ reads 

G^! m = J2 l.i W & i) ■ E^ - n« • n«>) (e£> ■ E% - n<£ ■ n»>)#(J, 0- (2) 

This expression has to be calculated using the double grasping operators acting on the nodes 
of the 4-symplex spinnetwork derived in |8] and a suitable choice of boundary state that we 
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want picked on the classical geometry q of a regular 4 simplex. The boundary state was 
defined [5] [8] [21] as gaussian wave packet in spin and intertwiners variables, centered on 
the values jo and io determined respectively by the areas and the angles of the background 
4-simplex. As in ordinary quantum mechanics a crucial role was played by the phase of the 
packet. The spin phase coefficients were fixed by the background extrinsic geometry [5]- The 
intertwiner phase coefficients were fixed by the requirement that the state remain peaked 
after a change of pairing to the value io determined by the classical value of the 3d dihedral 
angles of each regular tetrahedron inside the 4-simplex. The crucial point is that the non 
commutativity 23 of the the different angles of a tetrahedron, represented by the intertwiner 
variables in different pairings, requires a state with a phase dependence in the intertwiner 
variables to be peaked on the background angles in any pairing. The correct value |25l [8] for 
this in the equilateral case is exp{i^i n }. In [5] and [8], (a suitable adjustment of) the BC 
vertex was chosen for W and in this limit the propagator depends only on its asymptotic 
behavior, this has the structure [24] WbcQ) ~ e i (ijG ' 5j) e i *" 5j + e-iWe-'*' 41 where G 
is the 10 x 10 matrix given by the second derivatives of the 4d Regge action around the 
symmetric state, and $ is a lOd vector with all equal components, which were shown [5] to 
precisely match those determined by the background extrinsic curvature. Computing ^ with 
these ingredients the crucial point is that the phase in the link variable in the boundary state 
cancels with the phase of one of the two terms of Wbc, while the other term is suppressed 
[5] for large jo (this was the key mechanism that gave the correct diagonal components with 
the BC vertex) but the rapidly oscillating factor in the intertwiners variables is completely 
uncompensed by the dynamics and suppress the integral^. The intertwiner independence of 
the BC vertex prevents the propagator to have the correct long distance behavior. In Ref. 
[21] has been proposed an asymptotic form of a vertex W that includes a gaussian in all the 
15 variables, and most crucially a phase dependence also on the intertwiner variables. The 
proposed form for the asymptotic of W was W(j,i) = e t(«G«) e ^« + e -|(MGM) e -i*-H ^ 
where G is a 15 x 15 real matrix that scales as 1/jo- The quantity (j> — (<f> nm , (fin) is now a 
15d vector: its 10 spin components (j> nm just reproduce the spin phase dependence of the BC 
vertex; while its five intertwiner components are equal and fixed to the value <p n = ^. This 
phase dependence is the crucial detail, that makes the calculation work because it allows the 
cancellation of the phases between the propagation kernel and the boundary state through 
which the dynamical kernel reproduces the semiclassical dynamics in quantum mechanics. 
If this does not happens, the rapidly oscillating phases suppress the amplitude. The results 
of Ref. [8] motivated the search for alternative models. Ref. [21] shows that it is possible 
to recover the full propagator of the linearized theory from the LQG propagator and gives 
indications on the asympthotic behavior of an alternative vertex: in particular it requires for 
the new models an oscillation in the intertwiners. A natural question was then if the new 
models present or not the desired oscillation. 

Let us focus on the EPRL model [10]. For given Immirzi parameter 7, the vertex ampli- 
tude is defined as a function of five SO(3) intertwiners i n and ten spins j nm 

w( jnm ,i n ) = {i5jM |1 ~^'" m ,^) {i5iM (1+ ^" ro ,ffi Ufths^ • (3) 

The functions {15j}jv are normalized 15j-symbols and the are fusion coefficients from 

SO(3) to SU(2) L x SU(2) R introduced in [9]. Indeed the model differs from the BC one 
only for the structure of these coefficients. One of the main open question on these new 
models was their large large spin's asymptotics to repeat the propagator calculations, with 
particular focus on the new intertwiner dependence, and their semiclassical limit to establish 
a link with Regge calculus. The analysis of these models started numerically using a new 
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technique to test the spinfoam dynamics introduced in Q3]: The propagation of semiclassical 
wavepackets. This technique is based on the following simple consideration: In ordinary 
quantum mechanics the propagation kernel Wt(x,y) = {x\e~n Ht \y) of a one-dimensional 
nonrelativistic quantum system defined by a hamiltonian operator H can be used to check 
the semiclassical behaviour of the theory using the propagation of wave packets 4ix, P (x) with 
Wt(x,y). If the propagation is coherent, given a semiclassical wave packet tp Xi , Pi (y) centered 
on the initial values Xi,Pi its evolution under the kernel <j)(x) ~ J dy Wt{x,y) ^pxi, Pi (y) 
will be a new semiclassical wave packet centered around the final data Xf,pf. Now the 
equations of motion of any dynamical system can be expressed as constraints on the set 
formed by the initial, final variables namely on the boundary variables. This means that 
a solution of Einstein equations can be recovered by an appropriate set of constraints on 
the boundary variables. The derivation of the vertex amplitude presented in [9] indicated 
that the process described by one vertex can be seen as the dynamics of a single cell in a 
Regge triangulation of general relativity. It follows that a boundary wave packet centered 
on the boundary variables of a flat 4-symplex must be correctly propagated by the vertex 
amplitude evolving in a wave packet centered on the correct boundary variables, if the vertex 
amplitude is to give the Einstein equations in the classical limit. The tested process was the 
evolution of four identical coherent tetrahedra [25] with the propagation kernel of the EPR 
model. A single regular Regge cell has all the boundary variables equals: this means that if 
the quantum kernel describes the same process, the final state should have been a coherent 
tetrahedron with the same geometrical properties of the evolved states. In |14| the expected 
mean and phase were found and in [TS] were found the correct analytical properties of the 
evolved state. To calculate this propagation were used the asymptotic properties of the fusion 
coefficient [16] that allowed drastic simplifications in the numerical calculations and a new 
numerical algorithm (using techniques similar to [26]). The asymptotic properties of fusion 
coefficients showed that these coefficients not only give nontrivial dynamics to intertwiners 
at the quantum level, but also map semiclassical SO(3) tetrahedra into semiclassical SO(4) 
tetrahedra. To complete the analysis of the asymptotic properties of the model with this 
approach it is needed the asymptotic formula for the 15j-symbol that is still missing. Though 
was used the drastic approximation of fixing all spins, these were the first indication that the 
EPR model has the good semiclassical limit and the propagation scheme gave an hint on the 
required linear phase dependence in the asymptotic regime (in fact the phase of the outgoing 
state was exactly -| and this phase comes directly from the vertex) . In the meanwhile dealing 
with the tentative of relating the spinfoam models to a physical scalar product a new integral 
formulation has been given [27] . 

Key developments appeared in [18] were the authors studied the semiclassical properties 
of the FK model using the construction of a path integral with a discrete and local action 
introduced in 17 • They showed that, in the semi-classical limit, the amplitude converges 
rapidly towards the exponential of i Sn eg ge, if the faces spins can be understood as coming 
from a discrete geometry, otherwise the spin foam amplitude is exponentially suppressed. 
Remarkably this result holds for an arbitrary triangulation and gave important informations 
on the good semiclassical limit of the model and its relation with Regge calculus, but it 
was not yet viable for the explicit calculation of the graviton propagator, where the explicit 
dependence on the intertwiners had to be understood. 

A new crucial result is the semiclassical limit recently performed by Barrett, Fairbairn and 
collaborators [2D] and by Conrady, Freidel [TpJ that gives the exponential of the Regge action 
for the vertex projected on the Livine-Speziale coherent state basis. This important result 
is now directly applicable in the propagator calculations. In [28] it is noted that this result 
implies for the equilateral configuration WeprUo, io + Si) ~ Ne~ lS ^ ^ e 1 "? 1 . The important 
point here is the appearance of the correct ~ factor in the phase, which was missing in the 
BC vertex. Thus, the new vertex has the asymptotic behavior that was guessed in [2T], in 
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order to yield the graviton propagator. The obstacle that prevented the BC vertex to yield 
the proper propagator is resolved by the new vertex. The complete calculation with the new 
vertex is now possible [29] and in course. 
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